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Abstract 

For the first time we present analytical results for the contribution of electromagnetic fluctuations 
into thermodynamic properties of modulated systems, like cholesteric or smectic liquid crystalline 
films. In the case of small dielectric anisotropy we have derived explicit analytical expressions for 
the chemical potential of such systems. Two limiting cases were specifically considered: (i) the Van 
der Waals (VdW) limit, i.e., in the case when the retardation of the electromagnetic interactions can 
be neglected; and (ii) the Casimir limit, i.e. when the effects of retardation becomes considerable. 
It was shown that in the Casimir limit, the film chemical potential oscillates with the thickness 
of the film. This non-monotonic dependence of the chemical potential on the film thickness can 
lead to step-wise wetting phenomena, surface anchoring reorientation and other important effects. 
Applications of the results may concern the various systems in soft matter or condensed matter 
physics with multilayer or modulated structures. 

PACS numbers: 68.15. +e, 34.20.-b, 78.67.-n 
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I. INTRODUCTION 



Stratified systems are omnipresent in soft matter and sofid state physics. Understand- 
ing of molecular interactions in these systems is an important step for controlling their 
properties and aggregation processes. Though these interactions may have many differ- 
ent specific sources, the electromagnetic fiuctuations (VdW and Casimir forces) are their 
common cause. These interactions become appreciable in the submicron range and rapidly 
increase in the nanometer scale. Although, evidently, the interactions would not scale nor 
ad infinitum, neither at shorter scales, where steric short range interactions, i.e. interaction 
forces whose origin is in the Pauli principle, become dominating ones. Still, submicron - 
nanometer range, where long range forces are dominant ones, is the most important range 
for applications in nanotechnology processing. Several approaches have been employed to 
calculate these fluctuational forces. The simplest one is Hamaker approach [1], summing 
all pairwise interactions, which is evidently not applicable to condensed media, where all 
molecular units are strongly correlated. The rigorous continuum method was developed in 
[2] (see also [3]), by Lifshitz, Pitaevskii, and one of the authors of the present paper (I.D.) 
who derived the general expression for the electromagnetic fiuctuational interaction between 
two macroscopic isotropic bodies separated by an isotropic film. Although the generalization 
of the analysis [2] for the anisotropic and multilayer systems is conceptually straightforward, 
it deserves some attention, as it implies prohibitively tedious and bulky calculations, which 
could be done analytically only under certain rather restrictive approximations (see e.g., [4], 
[5]). Luckily, however, there are some cases when these approximative calculations can be 
performed in a well controlled way, and besides they may be useful for systems interest- 
ing for applications. For example when dielectric anisotropy is small in the whole region 
of frequencies relevant for the interactions, one can use the regular perturbation theory [6], 
which is a typical case for liquid crystalline films [7], [8]. Moreover, cholesteric liquid crystals 
having a long wavelength helical orientational modulation can be treated as a continuous 
model of multilayer systems. 

In what follows we calculate the electromagnetic fiuctuational contributions to the chem- 
ical potential of a cholesteric or smectic liquid crystal film (thickness I) confined between 
two semi infinite isotropic media (having in mind experimentally interesting cases, the lat- 
ter can be, e.g., solid, glass- like substrate and air). Our approach is a macroscopical one. 
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i.e. we consider quantities averaged over physically infinitesimal volumes. Thus we restrict 
ourselves to only the long wavelength part of the electromagnetic field fluctuations. In this 
range everything can be expressed in terms of macroscopic characteristics of the system 
(dielectric and magnetic permeabilities). 

Since for cholesteric liquid crystals typical scales of orientational modulations are very 
large (of the order of 500nm and even larger for nematic - cholesteric mixtures [7], [8]), 
they may be the perfect candidates for observing of predicted (by our macroscopic con- 
tinuum approach) behavior. With a certain modification in order to take into account 
magnetic properties, one can also perform analogous calculations for large pitch magnetic 
spiral structures. The same is true for so-called lyotropic smectic liquid crystals where den- 
sity modulations occur at scales larger than characteristic molecular lengths (2 — 5)nm), 
however, for thermotropic smectic liquid crystals this window of validity of our approach is 
much more narrow (if it exists at all). In this case the continuum theory loses all pretense 
of quantitative predictability and a complete description should incorporate microscopic 
structure and molecular short range forces what can be done only numerically (see e.g., [9], 
[10]). 

Our motivation for presenting this discussion is two new predictions which have emanated 
from our investigation. Namely, in the Casimir limit we have found spatially oscillating terms 
in the chemical potential of the film. In the VdW limit (small film thickness) we have found 
a monotonic in I finite anisotropic correction to the Hamaker constant, making possible a 
reorientation anchoring phase transition with variation of the film thickness. 

The remainder of our paper has the following structure. Section II contains basic me- 
thodical details and equations necessary for our investigation in the frame work of a regular 
method for calculating higher order perturbative corrections. New results are discussed in 
section III in subsections III A and IIIB, where we present the electromagnetic fluctuational 
contributions into the cholesteric film chemical potential. Finally, Section IV deals with 
miscellaneous subjects related to the physical consequences of our results for wetting, sta- 
bility and anchoring phenomena in cholesteric films. Some technical material is collected in 
Appendix to the paper, and those readers who are not very interested in derivations can 
skip this Appendix finding all essential physical results in the main text of the paper. 
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II. ELECTROMAGNETIC FLUCTUATIONS IN CHOLESTERIC FILMS 



We limit our analysis to the long range forces induced by electromagnetic fluctuations, 
thus neglecting all kinds of structural thermal fluctuations. We consider a leftmost half 
space, — oo < x < 0, as an isotropic body (denoted in what follows by the indices 1 (solid, 
e.g., glass substrate) and a rightmost half space 2, I < x < +oo, (vacuum or air), separated 
by non-homogeneous and anisotropic film, < x < I, (denoted by the index 3) with thickness 
I. 

First closely following the general method [2] we study the case when the substance 3 is 
a cholesteric liquid crystal film. We chose the system for our benchmark to determine all 
specific features of electromagnetic fiuctuational forces in anisotropic modulated systems. 
We take the x-axis perpendicular to the separating surfaces (planes x = and x = I). This 
is not the place to explain the general approach to fiuctuational electromagnetic forces in 
detailed, however, in a stripped down version the theory [2] is reduced to the calculation of 
Green functions Dik to the Maxwell equations, i.e., 

[eii{r,i\un\)ul- curli^curl^i]Dik{r,r';ujn) = -47r5(r - r')^^^ , (2.1) 

where ej;(r, i|u;„|) is dielectric permeability at Matsubara imaginary frequencies ic<j„. It is 
worth noting that the dielectric function is always real on the imaginary axis irrespective of 
whether the system is absorbing or not. In what follows we use atomic units, i.e. put light 
velocity c and Planck constant ^ as 1 (except where explicitely stated to the contrary and 
the occurrences of ^ or c are necessary for understanding). 

Calculating somehow D^k one can find the stress tensor in the film 

(^ik^-^J^Vi^rl ^mfc(r,r;ia;„) - -5mfc^g(r,r;ia;„) - -6ikDfj{r,r;iuj,,)[ . (2.2) 

n=0 ^ '- ' 

Here, a^^ is electromagnetic fiuctuation part of a^k, and functions Dff^ and Dfl are corre- 
sponding averages of fiuctuating electric and magnetic field components 

Df,^{r, r'; iujn) = -ujlDu{r, r'; iwn) ; D!^{r, r'; iujn) = curlucurl'^^Dimir, r'; iwn) . (2.3) 

Because the problem is homogeneous in z — y plane we make a Fourier transform with respect 
to z — z' and y — y' coordinates, thus getting Dii;{x,x' ,(i;iun). Taking into account that 
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4 = ei5ik, = 5ik, and 












63 + e„(l + cos(2goa:;)) 
Ea sin(2goa:;) 



ea sm(2goa:;) 
£3 + ea(l - cos{2qQx)) 



(2.4) 



where we define the cholesteric dielectric permeabihty as = ^z^ik + '^^a'^ii^k (ii — 
(0 , cos{qQx) , sin(go3:^)) is cholesteric equilibrium director, go is cholesteric spiral modula- 
tion wave vector, related to cholesteric pitch p ets qo — 27r/p, and in what follows we assume 
a small dielectric anisotropy < 63, what is the case for all known cholesteric hquid crystals 



In the natural coordinate frame x, y, z (2.2) is a set of nine coupled equations (presented 
in the appendix to the paper). This set should be supplemented by the usual boundary 
conditions (corresponding to the continuity of the tangential components of the electric and 
magnetic fields) which are the continuity conditions for 



In principle the full information concerning electromagnetic fluctuational forces is contained 
in solution of (Al) - (A9) supplemented by the boundary conditions (2.5). In practice, 
however, for a given cholesteric dielectric permeabihty (2.4) the analj^ic solution is infeasible, 
thus some approximations are definitely needed not only to make calculations possible, but 
also to understand correctly underlying physical phenomena. Luckily since < £3, for all 
known cholesteric liquid crystals, one can solve the equations using a regular perturbation 
theory with respect to the small parameter ea/e^. 

Omitting a large amount of tedious (although straitforward) algebra, we get in the first 
order approximation the set of four equations to find D^z , Dyy , D^y, and D^x- Two latter 
equations are simple relations 



where we designate = euj^ + q^, and as e one should take 1 , 62 and £3 + £„, for air, 
substrate, and cholesteric film respectively. Two other equations (for D^z and Dyy) deserve 
some precaution only in the cholesteric film region {0 < x < I) where they are equivalent 
to one particle Schrodinger equations in a periodic potential (see appendix, (AlO), (All)). 



[7], [8]. 



Dyk , Dzk , curlyiDik , curlziDik . 



(2.5) 




(2.6) 
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However, since to derive these equations we already used the perturbation theory, we have to 
consider only so-called almost free particle approximation (weak coupling limit) , which gives 
for our case (A12), (A13). Determining arbitrary constants Ai , Cj entering these expressions 
from the boundary conditions (2.5) and subtracting homogeneous system contribution, we 
find finally almost conventional [2] expressions for the Green functions 

An 



cosh[w3(a; — x')] , 



where, however, A contains oscillating terms 

A = 1 - exp(2^/;30 (^3 + ^i)(t^3 + ^^2) .-^ ^ ^.^^25^;) + cos(2go0] , 

(W3 - Wi){W3 - W2) 



and we designate 



and 



2{ql + wl) 



Q!2 — 



W3 _^ 



4(g2 + ^^2) ■ 



Analogously 



cosh['u;3(x — x')] , 



where 



Ai = 1 — exp(2'u;3Z) 



(w3 + W2)(eiW3 + Wi) 
{Ws - W2){eiW3 - Wl) 



and as above 



+ 



(^aW2W3qo 



A{qi + wl)es{w3 + W2) ' Ae^iqi + wD' 



and 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



[1 + /3isin(2go0 + /32cos(2go0] , (2-12) 



(2.13) 



/32 = 



(2.14) 



4e3(g^ + w|)(w3 + W2) 2e3{q^ + wl)' 
These expressions (2.7) - (2.14) are our main results and they are ready for further inspection. 
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III. FLUCTUATIONAL CONTRIBUTION TO THE CHEMICAL POTENTIAL 



Our aim is to calculate fluctuational contributions into thermodynamic properties of a 
thin cholesteric film on an isotropic solid substrate. Since the film should be in equilibrium 

P-Pch- , (3.1) 

where a'^^ is given in (2.2), P is the air or vapor pressure, and Pch is the pressure of the bulk 
cholesteric. Note that using a standard definition for the pressure 

Pch^P^-E, (3.2) 
dp 

were p is cholesteric mass density, and into the energy density E, orientational deformation 
energy (Frank energy [7], [8]) has to be added to isotropic liquid pressure Pq. Therefore, 

Pch = Pq + K22qol{eiknniV - Qo) , (3.3) 

where Eikn is antisymmetric tensor, K22 - twist orientational elastic modulus, and 7 = 
— {dlnqo/dp) is cholesteric pitch "compressibility". Now wc are in the position to reap 
fruits of our calculations. From (3.1) and given above expressions (2.2), (2.7) and (2.11) the 
electromagnetic fluctuational part of the cholesteric film chemical potential p{l) reads 

A'O-^E/ '<'«»3(^ + ^). (3.4) 

Putting all expressions together one can write down (3.4) explicitely. This general expression 
is a very bulky one, however, it can be considerably simplified in two important limiting 
cases. 



A. VdW interactions 

We term VdW interaction the case of small thicknesses, i.e. / <^ A*, where A* characterizes 
the main absorption band of the cholesteric film. Typically in known cholesteric liquid 
crystals [8], A* is of the order of 10~^cm and qoX* < 1. Thus we can perform an expansion 
over qqI in (2.7) - (3.4) and get a relatively compact answer for this limit 

p{l)^Po{l)+a^, (3.5) 
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where /io(0 ocl ^ is known [2] for isotropic films VdW contribution, and 



dp dC 

Jo Jo 



1 -1 



a — 



63(^0 Le3(^0-e2(^0)(e3(^0-l) 



exp(x) — 1 



(3.6) 



Let us emphasize in passing that this anisotropic VdW contribution scales as l/l, and it is 
more long ranged than known fluctuational forces in isotropic homogeneous materials. It 
comes from the interplay of the characteristic lengths (film thickness and cholesteric pitch) , 
and the extra factor in (3.5) provides correct dimensionality. Note also that the sign of 
this contribution (3.6) does not depend on the cholesteric spiral chirality (left or right) as 
one can certainly expect but it does depend on dielectric properties and could be positive 
or negative. In the latter case the planar director orientation (cholesteric axis along x- 
axis) is always stable, and in the former case this orientation is energetically unfavorable 
and could become unstable at a certain thickness l^r- The critical thickness is determined 
by a competition of VdW contribution (3.6) and stabilizing planar orientation short range 
anchoring energy W 



This transition is analogous to the Preedericksz transition known in liquid crystals [7], [8], 
however, while the Preedericksz transition in liquid crystals is driven by the quadratic cou- 
pling between an external magnetic (or electric) field and the director, in our case the 
transition is induced by internal electromagnetic field fiuctuations (VdW forces). 

Some comments about the validity of (3.7) seem in order here. The critical thickness (3.7) 
is determined by a competition of the VdW (long range) contribution into the bulk chem- 
ical potential, and phenomenological (basically short range) anchoring energy W. Clearly 
the approach misses short range contributions into the bulk properties. If the only forces 
between particles had a range of the order of molecular dimensions (short range forces) , the 
corresponding contributions would decrease exponentially with increasing distance, therefore 
assuming that at the scale (3.7) the VdW contribution is dominating, we believe that Icr is 
much larger than all molecular scales. This assumption is not always justifiable, and in this 
case, (3.7) does not maintain numerical accuracy, but nevertheless is useful and instructive 
for gaining some qualitative insight on interplay short and long range forces. 




(3.7) 
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B. Casimir force 



We now turn to the opposite limiting case I > X* when electromagnetic interaction 
retardation effects become relevant, and it is referred as Casimir limit. Closely following to 
[2] we can replace in this case all dielectric permeabilities entering the general formula (3.4) 
for by their static values, and after some algebra end up with 

lio{l) + b^cos{2qol), (3.8) 

where /^o oc is conventional isotropic Casimir contribution, and the coefficient b reads as 

1 foo f-oo a;3 0cxp(,T) , , 



where we designate 



and 



(gio+j^a(o))(,s2o + p) 
(sio -pei(0))(s2o -p) ' 



(3.10) 



.o = ^|||-l+P^.o = y^-l+P^ (3.11) 

The Casimir limit (3.8) is particularly instructive, because it shows oscillating with a thick- 
ness dependence of the chemical potential. Since the both contributions (monotonic and 
oscillating ones) decay as the same power of Z, it is especially interesting to hypothesize 
the case (not forbidden in principle although formally beyond our first order perturbation 
theory approximation) when it becomes negative. Suffice it to say that such an oscillating 
function //(/) leads to a ffim formation which is stable only in a certain range of its thick- 
ness. This kind of successive thinning or wetting - drying transitions are known for layered 
smectic phases (see e.g. [11], [12], [13]). We have shown that Casimir forces in cholesteric 
liquid crystals can be responsible for similar phenomena. Lacking sufficient data on dielec- 
tric permeability frequency dispersion we can at present discuss only the general qualitative 
features of Casimir forces in cholesteric liquid crystals. 

Actually our picture is not entirely correct since we have neglected all short range forces, 
and this relatively weak oscillating Casimir contribution can be swamped by stronger short 
range forces. However by their nature, the short range forces have no connection with the 
effects under consideration which are due to cholesteric long wavelength orientational modu- 
lations. Short range contributions to the chemical potential are the same in cholesterics and 
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in analogous homogeneous systems (e.g., nematics). Thus, even in the case of dominating 
short range forces, the Casimir contribution (3.8) could be disentangled in differential mea- 
surements. We anticipate that the discussed above effects of Casimir film instabihties will 

be observable and that understanding of underlying mechanisms will be essential to predict 
and to use these phenomena. 

IV. CONCLUSION 

In summary, in this paper we have calculated electromagnetic fluctuational contributions 
into cholesteric liquid crystal film chemical potential, and found the oscillating with the 
film thickness contribution (3.8). Quite remarkably this result illustrates that the collective 
nature of the Casimir forces suggests it to have a non-trivial dependence on film thickness, 
and even the sign of the Casimir force is dependent of material parameters and structure. 

Our theoretical approach can be extended in several directions. First, our main result 
(3.8) with evident modifications can be applied to a smectic liquid crystal film, where its 
dielectric permeability acquires an oscillating contribution owing to one dimensional density 
modulation. Indeed, for a smectic film we should replace (2.4) by 



where we chose x-axis as a normal to smectic layers, and Qsm is smectic density modulation 
wave vector. Of course the macroscopic approach, we used in this paper, is valid only when 
the corresponding modulation periods are large with respect to atomic or molecular scales. 
It is always the case for cholesteric liquid crystals (a typical pitch is about (500 — 700)nm), 
but usually not true for standard thermotropic smectics, where their density modulated 
with a period in the range of (2 — 5)nm [8], [7]. Thus as an application of our results 
to smectics, we have in mind long period lyotropic smectic liquid crystals and some other 
lamellar (layered) membrane structures [14] where period can be much larger. 

One more comment of caution. The relation between the fiuctuational force and chemical 
potential (in fact the mechanical equilibrium condition of the type of (2.6) we have used for 
cholesteric films), is valid for liquids (or uniform in density cholesteric liquid crystals having 
no solid-like elasticity). In smectics, like in solids, an elastic layer compressibility modifies 
the equilibrium, and there is no simple relation between electromagnetic fiuctuation force 




(4.1) 
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and film chemical potential. For these cases our theoretical predictions should be confronted 
with direct force or (in anisotropic cases) angular torque measurements. 

On equal footing one can consider an apparently unrelated problem of electromagnetic 
fluctuations in solid long periodic magnetic structures. The issue of Casimir forces in mate- 
rials with nontrivial magnetic susceptibility has been discussed recently [15], [16], [17]. To 
treat the magnetic case we have to include magnetic permeabilities /i of all three media into 
the equation for the Green function (2.1) and to modify the boundary conditions (2.5). The 
latter ones are read as the continuity of 



^ZZ -I 



1 dD, 



) ^yy 1 



e dD. 



yy 



(4.2) 



jji dx ' ^ dx ^ 
and besides w functions are now -u;^ = e^cj^^+g^. Calculations become more involved (see e.g. 
[18]) but in the frame work of perturbation theory are still feasible analytically. However 
results and even signs of fluctuational contributions depend on many unknown functions 
describing dielectric and magnetic permeability frequency dispersions. A more specific study 
might become appropriate should suitable experimental results become available. 
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APPENDIX A 



We take the y-axis along the vector q and in the coordinate frame attached to the system 
(2.2) in self-evident notations reads 



{^zzUJn + 



dx^ 



Dzz + ^zy^n^yz = -^'k5{x - x') , 



(Al) 



^2 



d 

Dyy + ey^LolD^y ^ i—qD^y = -4:7id{x - x') , 



(A2) 



d 

{^xxi^l + (f')Dxx + iQ^Dy^ = -47r5{x - x') , 



dx 



dx^ 



Dzy + ^zy^n-^yy — ) 



(A3) 
(A4) 
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0, 



(A5) 



2 9 
Dyx + ^eyzU^D^x + iq—D^j. 



0, 



(A6) 



d 2 

Dyz + iq^D^z + e^yU^D^z = , 



(A7) 
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(A8) 
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{^xx^jI + q^)Dxy + iq-^Dyy = . 



(A9) 



From these equations in the first order over the small parameter ea/ea we get two decoupled 
equations for D^z and D^x, which are in the only non-trivial region < x < I 



+ eaUJn cos(2goa;) 



Dzz = —4:7r6{x — x ) , 



(AlO) 



and 



w 



2 9^ taWl taq^ 



^yy ~ ) 



(All) 



2 dx^ 63 

and from the solution to these equations using two relations (2.6) one can find two others 



Green functions Dxy and Dxx- General solutions to the equations (AlO), (All) read as 



Dzz = Ci exp{w3x) 



4goW3 



sin(2goic) 



4w3(w| + q^) 



{ws cos{2qox) + qo sin(2goa:)) 



+C2exp{—W3x) 



1 + sin(2goa^) - rY^, — 2t(~^3 cos(2go2^) + qoSm{2qox)) 



27r , 

exp(— kr — X 

W3 



4goW3 

1 + 



^qo{wl + ql) 



4w3(w| + ql) 
qo cos{2qox) + W3 sin(2goa:;) 



(A12) 



and 

Dyy = Ai exp{w3x) 



1 + sin(2goa;) - - — t^-^ — 2\ (""^3 cos(2goa;) + qo sin(2goa;)) + 



4goe3 



+A2exp{-W3x) 



1 - sm(2goa:;) + 

463^0 



463(^0+^1) 



4^3(^1 + qi) 



W3 cos{2qox) + qo sin(25oa;)) + 



27rw3 , TTW^ea , . 

^ exp(-W3|x - X \) exp[w3X)t+ 



e3< 



H — exp(i(;3x)s+ + ^ ^ exp(-W3a;)p+ - ^. .9 exp(-W3a;)m+ , 



3"^n 



3'^n 



463 
463^1. 



(A13) 
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where Cj , Ai are constants which should be found from the x — and x — I boundary 
conditions, and in (A13) we introduce shorthand notations 

, -W3cos{2qox) + qosm{2qox) 1 / / o ^ 

t+ = exp{w3X - 2wzx) ^ — 2^ ; s+ = -- — ex^{W2,x - 2w2,) ; 

1 1 

m+ = exp(— wsx' + 2w^x) ; p+ = exp(w3a;') — sin(2goa:^) • (^14) 

2^3 2gfo 
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